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Abstract 

In the framework of the teleparaUel equivalent of general relativity the energy 
density Eg of asymptoticaly flat graviational fields can be naturally and unambigu- 
ously defined. Upon integration of Eg over the whole three dimensional space we 
obtain the ADM energy. We use £g to calculate the energy inside a Schawrzschild 
black hole. 
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The definition of localized energy density is a longstanding problem in the theory of general 
relativity |T]]. The variation of the Lagrangian density with respect to the metric field yields 
the energy-momentum tensor of the standard field theories. It is well known that this 
procedure cannot be satisfactorily applied to the Hilbert-Einstein action. Usually it is 
asserted on the basis of the principle of equivalence that the gravitational energy cannot 
be localized^]. The principle of equivalence is frequently invoked to assure that the 
gravitational field can be made to vanish in a sufficiently small region of the spacetime: 
the dynamics of a system under the action of a gravitational field in a locally inertial 
frame remains unchanged if we eliminate the gravitational field but consider instead an 
appropriate non-inertial frame. For this reason there have been attempts to define quasi- 
local energy in general relativity (see ref. and references therein). 

The status of the principle of equivalence in the formulation of general relativity, as 
a basic principle of the theory, is, to some extent, controversial (see the discussion in the 
Preface of Synge's book^]). Einstein's principle of equivalence^] estabhshes an equiva- 
lence between inertial and non-inertial reference frames, from what follows the equality 
between inertial and gravitational masses, and ultimately the prescription according to 
which the gravitational field couples in the same way with the other fields and matter 
systems in nature (considering, of course, the different couplings between gravity and 
boson and fermion fields). However, the symmetry group in general relativity is the 
group of transformation of coordinates, and a coordinate transformation can have no ef- 
fect on the presence or absence of a gravitational field. It is precisely the requirement 
of invariance under general coordinate transformations that leads to the Einstein tensor. 
Cartan[^] proved that the most general second rank tensor (i) constructed in a coordi- 
nate independent way from the metric tensor and its first and second partial derivatives, 
{ii) having a vanishing divergence, and [Hi) linear in the second derivatives is Einstein 
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tensor (Lovelock showed that one can dispense with {in)). From the mathematical 
point of view Einstein's principle of equivalence plays no role in the determination of the 
gravitational field equations in vacuo. 

Therefore the principle of equivalence cannot be taken to rule out the existence of an 
expression for the gravitational energy density. The absence of the latter in the literature 
is really due to the description of the gravitational field in terms of the Hilbert-Einstein 
action integral, which is not an appropriate framework for such considerations. 

In this paper we will show that in the context of the teleparallel equivalent of general 
relativity (TEGR) there is a natural and simple definition of energy density of the grav- 
itational field for asymptoticaly fiat spacetimes. This expression is first noticed in the 
Hamiltonian formulation of the TEGR; specifically, it arises in the integral form of the 
Hamiltonian constraint. However, it also appears as a surface term in the total Hamilto- 
nian of the TEGR, and therefore it has a canonical significance. 

Notation: spacetime indices fi,^, ... and local Lorentz indices a,b, ... run from to 3. 
In the 3-1-1 decomposition latin indices from the middle of the alphabet indicate space 
indices, according to fi = 0,i, a = (0), (z). The tetrad field e'^ ^ and the (arbitrary) 
spin affine connection u^ab yield the usual definitions of the torsion and curvature tensors: 
R'^bf^u = d^Uy"'b + cUJu^'b - T"- = d^fi"" y -hcj^^bC^- .... The fiat spacetime 
metric is fixed by ?7(o)(o) = — 1- 

Let us briefly recall the TEGR|^. The Lagrangian density in empty spacetime is given 



by 



L(e, a;. A) 



keT^^^^'Tabc + eX^'^^Rab^.Auj) • 



(1) 



The constant k equals 
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where G is the gravitational constant; e 



deti^e"" ^) and 



167rG' 



{A"*'^'^} are Lagrange multipliers. The tensor E"'"^ is defined as 
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•^abc ^ ^rpabc _|_ rpbac r-pcah\^ _|_ ^ (^fj^^J^^ ff^T^^ 
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which yields 



1 1 

Sabcrp rpabcrp I rpabcrp rparp 

-Labc — -l-abc + -^-L -l-bac — -L -L a , 

where = T^q. 

The affine connection u^ab can be (identically) rewritten as u^ab = "i^^iab + K^ab, 
where "oj^afe is the Levi-Civita connection and A'^afe = ^^a ^^b "{Tx^j^^ + Tj,a^ — T^y\) is the 
contortion tensor. Such decomposition of uj^ab allows us to obtain the identity 

eR{e,uj) = eRie) + eS^^'^^T,^, - 2<9^(er^) , (2) 

by just substituting Uf^ab for °uj^ab + -f^/^afe on the left hand side (LHS) of (2). We observe 
that the vanishing of implies the equivalence of the scalar curvature R{e) - 

constructed out of e" ^ only - and the quadratic combination of the torsion tensor (for 
asymptoticaly fiat spacetimes the divergence in (2) does not contribute to the action 
integral) . 

Except for the presence of —k in (1), the Lagrangian density considered here is the 
same as in ref. [Q. Variation of L with respect to X"-'''^'^ yields 

R\^.uH = . (3) 

Let = denote the field equation satisfied by e°^. With the help of (3) it can be 
shown by explicit calculations that 



l{Ra,{e) - lea,Rie)} , (4) 



3 



again by replacing u^ab for °a;^afe(e) + K^ab in the LHS of (4). Therefore e"^ does satisfy 
Einstein's equations. 

The field equations arising from variations of L with respect to ujfj_ab can be best 
analysed in the Hamiltonian formulation. The latter has been presented in ref.|P], with 
the gauge cuoafe = being fixed from the outset. In this paper we will likewise maintain 
this gauge fixing, as it can be shown that in this restricted context the constraints of 
the theory constitute a first class set. The condition uoab = is fixed by breaking the 
local Lorentz symmetry of (1). We still make use of the residual time independent gauge 
symmetry to fix the usual time gauge condition e(fc) ° = e(o)j = 0. 

The Hamiltonian density H can be successfully constructed out of (1) in terms of 
canonical field variables and Lagrange multipliers. Because of the gauge fixing uoab = 0, 
H does not depend on p^"-''^ the momentum canonically conjugated to ujkab- Thus arbitrary 
variations of L = pq — H with respect to P'^"^ yields Ukab = (see expression (21) of the 
Lagrangian density in [8]; had we not fixed ujoab = the corresponding equation would 
be Ukab — DkUjQab = 0, which is equivalent to the Lagrangian field equation Rabok = 0). 
Therefore in view of uj^ab = 0, Ukab drops out from our considerations. The above gauge 
fixing can be understood as the fixation of a global reference frame. 

The details of the 3+1 decomposition of (1) are given in [8], except that the surface 
terms in eq. (6) below were not considered in the latter. The canonical action integral 
becomes 




(5) 



H 



NC + N'C, + E^„n 



(6) 
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and A^* are the lapse and shift functions, and S^n = ~'^nm are Lagrange multiphers. 
The constraints are defined by 

C = dj{2keT^) - keJ:''^Tkij - ' l^') , (7) 

Ck = -e(,)fca,n(^')^ - n(^>r(,),, , (8) 

with e = det{e(^j)k and T* = g'^^e'^^^^T(^j)ik- We remark that (5) is invariant under global 
SO (3) and general coordinate transformations. An important feature of this framework is 
that although we are considering asymptoticaly flat gravitational fields, the action integral 
determined by (1) does not require any additional surface term, as it is invariant under 
coordinate transformations that preserve the asymptotic structure of the field quantities. 
A clear discussion concerning the necessity of the addition of a surface term to the Hilbert- 
Einstein action AhE) in the case of asymptotically fiat gravitational fields, is given by 
Faddeev0. 

We consider now that for r — * oo we have e(j)fc ~ r]jk + |^ifc(^)- In view of the relation 

/ d?xdj{eT^) = — ^ / dSkidihk - dkhii) = Eadm , (9) 
87rG J IottG Js 

where the surface integral is evaluated for r — >■ oo, we observe that the integral form of 
the Hamiltonian constraint C = may be rewritten as 

J d''x^^ke^:''m^J + ^(n^'n,^ - ^n^)} = ^adm ; (lo) 

the integration is over the whole tree dimensional space. Since dj{eT^) is a scalar density, 
from (9) and (10) we are naturally led to define the gravitational energy enclosed by a 
volume V of the space as 
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The expression above is manifestly invariant under general coordinate transformations of 
the three dimensional space-like hypersurface E and yields Eadm if the integration is over 
the whole E. 

Let us note from (6) that / (PxH evaluated from a set (e(j)A;, 11*^^)') that satisfy the 
field equations, in a coordinate system such that for r — oo we have N = 1, Nj = 0, also 
yields Eadm- 

The appearance of a scalar and vector densities in (6) is intimately related to the fact 
that it is not necessary to add a non-covariant surface term to (1). The surface term 
that must be added to the Hilbert-Einstein action is not well behaved under coordinate 
transformations, but as Faddeev stressed, the action integral is invariant]^]. Since it is 
precisely this surface term that leads to Eadm, it is not possible to obtain a satisfactory 
energy density out of Ahe- 

A similar difficulty occurs in the context of the Einstein-Cartan formulation. The 
Lagrangian density for asymptotically fiat gravitational fields in the first order differential 
formulation is given by 

Lec = ee''^'e"'Rabf.,{uj) - 2d^{ee^^e"'u,ab) ■ (12) 
After performing a canonical 3+1 decomposition of Lec, the Hamiltonian density acquires 



a surface term which reads 10 



e,^ = d,i2Nee^'e''^°iOjab - iV^e"*n„,) . (12) 
In the appropriate limit / d^xSj^^ yields Eadm, but no energy density can be defined from 
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The gravitational energy-momentum four- vector of theories with local gauge symmetry 



has already been discussed in the literature, in the context of Poincare gauge theories [pT 
Such analyses are always carried out in the Lagrangian framework, and can equally well 
be applied to the TEGR. We note, however, that the present considerations are derived 
from the Hamiltonian formulation. To our knowledge, the appearance of the ADM energy 
in the Hamiltonian constraint has not been previously noticed. 

Brown and YorkQ] have recently furnished an expression of quasilocal gravitational 
energy density for compact geometries, by resorting to a Hamilton- Jacobi-type analysis. 
They consider the action integral evaluated from field quantities that solve the classical 
equations of motion, and as a function of the time interval. Then the action A satisfies 
the Hamilton- Jacobi equation E = — where E expresses the energy of the classical 
solution. In analogy with this procedure they define the quasilocal energy associated with 
a space-like hypersurface S as minus the variation of the action with respect to a unit 
increase in proper time separation between the boundary 5 of E and its neighboring two- 
surface, as measured orthogonally to S at B. In the present case we can similarly define 
the energy density as minus the variation of the action with respect to the proper time 
N{x). For a given set of solutions of the classical equations of motion the energy density 
Eg can be defined as 

in agreement with (11). 

The fixation of the time gauge condition prevents the construction of a global S0(3,l) 
energy-momentum four vector P°^, because the space and time components of P° can no 
longer be mixed. Nevertheless we can identify the components of such a vector. The 



field quantities {n"*^} are defined by H"*^ — t^- H^^)^ is the true momentum canonically 
conjugated to e{j)h- However, because of the time gauge condition, in the 3+1 decompo- 
sition n^^'"^ must be expressed in terms of canonical field variables. One ultimately finds 
n(o)fc ^ 2eT'= (see the discussion following eq. (19) of [8]). 

We assume now that e(j)fe and II^-'^'^ satisfy the constraints. Prom Cj = we obtain 

5.nM^ = -e^'^^^Xl^^'Tij^ik ■ (15) 
Therefore we define the energy-momentum quadruplet P" as 



P« = 4n'^'= . (16) 

Since the constraints are satisfied we have 



pirn) _ _^_gMfc]-]-(i)i J. 



(18) 



In the following we will specialize the Hamiltonian formulation to the spherically 
symmetric geometry, in order to compute the gravitational energy inside a Schwarzschild 
black hole. We fix the triads e(^k)i as 



/ sin9cos(f) r cos9cos(j) —r sinOsincf) \ 



(19) 



sinOsincj) r cosdsincj) r sinOcoscj) 

\ cos9 —r sinO / 

where A = A(r, i); {k) is the line index and i is the column index. The one form e^'^) is 

defined as 
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from what follows 



= e^*^) rdr + e^^^ edO + e^^^ ^d(f) , 



e('=)e(fc) = e^^dr^ + dO^ + r'^ sin^ 9 dcf)^ . 

We also obtain e = det{e(^k)i) = r'^ sinOe^. For r ^ oo we require A(r) ~ 0{^). 

The symmetry reduction is performed directly in the Hamiltonian. We first determine 
the Killing vectors ^ of gij — e^^^ i^{k)3- Next we require the vanishing of the Lie derivative 

Li_{e~^Iiij) = 0, where liij = gjke(^ijill'^^^'^. We obtain 

e-^Uij = diag{A{r,t),B{r,t),B{r,t)sen^e) ] (20) 

A{r, t) and B{r, t) are arbitrary functions. Prom (20) we can calculate all {n^*^)^}. Upon 
substitution of the latter and (19) into (5) we find out, as expected, that there is no 
canonical field quantity conjugated to B{r, t). Thus we enforce B{r, t) = 0, which implies 
in n^'^)^ = n^'^)^ = 0. Defining H by H = fcr^ A and integrating over angles we finally 
obtain the action integral 



A ^ An J dtdr{U\ - NC - N^Ci} , 



:2ii 



C = 2ke\l - - ^e-n^ + , (22) 



Ci = -e'|:(e"'n) , (23) 



d 

e = IGTr/c— [r(l - e"^)] , (24) 
or 



The constraints C2 = Cg and C3 = vanish identically. 

The Hamiltonian formulation established by (21-24) is completely equivalent to the 



corresponding construction in the framework of the ADM formulation ||T2[, as it can be 
shown that the Hamiltonian and vector constraints are equivalent in both cases. If we 
choose a coordinate system such that A^^ = 0, then the constraints and the evolution 
equations for A and 11 yield the Schwarzschild solution, 



e-^ = 1 _ ^ , (25) 
r 

together with A^^ = e~'^^. 

The total energy associated with (25) may be calculated from the surface term (24). 
In analogy with (9) we obtain, as expected. 



Eadm = J dreg{r) = m. (26) 
We can also compute exactly the gravitational energy inside a black hole. It is given by 

f2mG 

Ebh = / drsgir) = 2m, (27) 



which implies that the total energy exterior to the black hole is —m. 

It is of interest to evaluate, in addition, the the gravitational energy inside a spherical 

surface of arbitrary radius R: 

Eg = (28) 

This is exactly the expression found by Brown and YorkQ] in their analysis of quasilocal 
gravitational energy of the Schwarzschild solution. The method developed by Brown 
and York, however, does not seem to be applicable to an arbitrary metric field. Problems 
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appear in the calculation of the quasilocal energy in the framework of the Kerr metric 1 13 1. 
On the contrary, given the triad components restricted to a three dimensional hypersurface 
of the Kerr type we can easily calculate Eg by means of (14). 

In conclusion, we find out that in the TEGR there is a natural, consistent and unam- 
biguous definition of gravitational energy density. This fact indicates that this framework 
is suitable for the Hamiltonian analysis of the gravitational field. In particular, the in- 
tegral form of the Hamiltonian constraint, eq.(lO), may become an energy eigenvalue 
equation in the canonical quantization program. 
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